Genetic evolution of Waldenström macroglobulinemia across disease stages Free

Waldenström Macroglobulinemia (WM) is an indolent B-cell lymphoplasmacytic lymphoma characterized by heterogeneous clinical behavior and variable progression patterns. While MYD88^L265Pand CXCR4 mutations are frequently encountered, the broader genomic landscape underpinning disease evolution from smoldering WM (SWM) to active treatment naïve WM (TNWM), post-treatment remission WM (PTWM) and relapsed/refractory WM (RRWM) remains incompletely understood. To address this gap, we studied the mutational profile across the WM disease spectrum.

We performed whole exome sequencing (WES) on sorted (CD19 and/or CD138) bone marrow aspirates of patients across the spectrum of WM. Paired germline sequencing was performed using peripheral blood mononuclear cells (PBMCs). Unique molecular identifier-tagged reads were deduplicated using fgbio. WES data were processed through the Sarek pipeline. Any valid mutation from either caller, present on the index or in COSMIC was considered. Copy number variations (CNVs) were inferred via CNVkit and summarized at the gene and arm level (≥50% directional change). Mean tumor and germline coverages were 174X and 74X, respectively. Fisher's exact test was used for enrichment and hypergeometric test for overlap analyses.

A total of 78 high-quality WM samples [15 SWM, 44 TNWM, 13 RRWM, 6 PTWM in remission (all with paired TNWM samples)] were sequenced. The mutational burden varied across disease states (median SWM=1.48 m/Mb, TNWM= 1.12 m/Mb, RRWM = 1.59 m/Mb; p=0.04). Density plots revealed bimodal VAF distributions, consistent with clonal and subclonal architectures in WM.

Overall, 77% harbored a MYD88^L265P mutation by WES. Among MYD88^WT^ patients, AS-PCR identified MYD88^L265P in 4 (25%) not detected by WES, raising the overall mutation rate to 83%, plus one with a non-L265P variant. Low VAF (<10%) MYD88^L265P mutations were seen in 6% of matched PBMCs. MYD88 copy number gain was seen in 2 (17%) MYD88^WT TNWM patients. CXCR4 mutations were detected in 38% by WES (18% each frameshift and nonsense). Two (4%) CXCR4^WT samples were positive by PCR, giving an overall rate of 40%. CXCR4 mutations were more frequent in RRWM (54%) vs TNWM (41%) and SWM (27%), though not statistically significant. Del(6q) and Del(17p) were noted in 16% (n=12) and 7% (n=5) of patient samples. All del(17p) events were noted in TNWM samples and in one paired post-treatment sample. Overall rate of del(17p) and TP53 mutation was 14% (n=11).

RRWM was enriched in ARID1A alterations, a gene implicated in chromatin remodeling and plasma cell identity formation [38% (n=5) in RRWM vs. 11% (n=6) in TNWM and SWM], compared to TNWM and SWM p=0.03]. Alterations in NF-kB signaling gene MAP3K14 were more frequent in RRWM (31%, n=4) compared to the SWM/TNWM cohort [(7%), n=4; p=0.036). We noted an enrichment in the TP53 mutation rate in RRMW [n=3, 23%, OR 7.7] compared to TNWM and SWM cohorts [n=2, 4%), p=0.04]. Mutations in FUBP1 (p=0.03), an upstream regulator of MYC, and CARD11 (p=.043), a B-cell signaling scaffold protein were also more frequent in RRWM. Notably, none of the SWM patients harbored an ARID1A alteration (p=0.057), suggesting acquisition during progression. CNV analysis revealed that MAP3K14 amplifications were frequent in RRWM (n=3, 23%) compared to SWM/TNWM [n=2 (4%), p=0.043].

RRWM samples showed increased CD79B mutations (15%) and significant overlap with the C5/ABC diffuse large B-cell lymphoma (DLBCL) genomic cluster (Chapuy et al., Nat. Med. 2018; p=0.007). This overlap remained significant after excluding MYD88^L265P (p=0.019), driven by CARD11, CD79B, and IGLL5, suggesting that a subset of RRWM may acquire an ABC DLBCL-like mutational profile. These findings warrant further investigation using transcriptional profiling.

This study provides a genomic map of WM evolution from precursor to relapse. Increasing somatic mutations and CNVs support stepwise progression. Enrichment of genes involved in chromatin remodeling and plasma cell identity formation (ARID1A), NF-kB or BCR signaling (CARD11, MAP3K14), and DNA damage repair (TP53) were noted at relapse, all of which have been shown to drive aggressive disease biology. Emergence of an ABC-like genotype in a proportion of RRWM indicates convergent evolution with aggressive lymphomas. Transcriptional studies are warranted to explore this further.